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ABSTRACT 

Context. IC348 is a nearby (~ 310 pc), young (~ 2 — 3Myr) open cluster with > 300 members identified from optical and in- 
frared observations. It comprises Young Stellar Objects in various evolutionary phases from protostars over disk-bearing to diskless 
pre-main sequence stars. This gives us the opportunity to study evolutionary effects in the high-energy emissions in a homogeneous 
environment. 

Aims. We study the properties of the coronae of the young low-mass stars in IC 348 combining X-ray and optical/infrared data. In 
particular, we intend to shed light on the dependence of X-ray luminosity and spectral hardness on evolutionary stage and on stellar 
parameters such as mass, effective temperature and bolometric luminosity. 

Methods. The four existing Chandra observations of IC 348 are merged, thus providing a deeper and spatially more complete X-ray 
view than previous X-ray studies of the cluster. We have compiled a comprehensive catalog of IC 348 members taking into account 
recent updates to the cluster census. Our data collection comprises fundamental stellar parameters, infrared excess indicating the 
presence of disks. Ha emission as a tracer of chromospheric emission or accretion and mass accretion rates. 

Results. We have detected 290 X-ray sources in four merged Chandra exposures, of which 187 are associated with known clus- 
ter members corresponding to a detection rate of ~ 60 % for the cluster members of IC 348 identified in optical/infrared studies. 
According to the most recent spectral classification of IC 348 members only four of the X-ray sources are brown dwarfs (spectral type 
M6 and later). The detection rate is highest for diskless Class III stars and increases with stellar mass. This may be explained with 
higher X-ray luminosities for higher mass and later evolutionary stage that is evident in the X-ray luminosity functions. In particular, 
we find that for the lowest examined masses (0.1 — 0.25 Mq) there is a difference between the X-ray luminosity functions of accret- 
ing and non-accreting stars (classified on the basis of their Hq emission strength) as well as those of disk-bearing and diskless stars 
(classified on the basis of the slope of the spectral energy distribution). These differences disappear for higher masses. This is related 
to our finding that the I/x/iboi ratio is non-constant across the mass/luminosity sequence of IC 348 with a decrease towards lower 
luminosity stars. Our analysis of an analogous stellar sample in the Orion Nebula Cluster suggests that the decline of Lx/iboi for 
young stars at the low-mass end of the stellar sequence is likely universal. 

Key words. X-rays: stars - stars: coronae, activity, pre-main sequence 
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1 . Introduction By means of optic al and IR spec troscopj y and ph otome- 

ti-y, the investigations oflLuhman et aP (ll998l) : lLuhmanl (119991) : 
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The open cluster IC348 is arguably the best investigated rich iLuhman et alll2U03l |2005a|) progressively increased the com- 
(> 100 members), and very young (~ 2 - 3 Myr old) stellar pietness of the census of the young stellar population in IC 348. 
cluster within - 300 pc from the Sun. IC 348 is associated to the These studies yielded a sample of 302 well characterized mem- 
Perseus molecular cloud complex at a distance of about 3 1 pc bers, for most of which individual stellar masses and ages can be 
dHerbig 1998) and located just to the south of the optically bright reliably estimated. Although the cluster contains a few embed- 
star o Per (Bl III). The cluster members spread over about 15' x jed stars this extinction-limited {Ay < 4 mag) sample of mem- 
15' on the sky. The optically brightest and most massive cluster bers is nearly complete down to masses of M > 0.03 Af©; it also 
member is the B5 V star BD H-3 1 643. contains 23 spectroscopically identified young brown dwarfs 
In a lai-ge number of observational studies, more than 300 in- (BDs), with spectral types as late as M9 and estimated masses 
dividual cluster members have been identified and well charac- down to 0.015 Mq. Recently, iBurgess et al.l ^OW^ discovered 
terized so far. The pioneering study of Herbig ( 1954) discovered a young T-dwaif candidate in IC 348. If its status as a cluster 
16 Ha emitting stars in IC 348. Later, Herbig ( 1998) identified member is confirmed, the derived spectral type of T6 suggests 
about 100 members by optical spectroscopy and photometry as a mass of only a few Jupiter masses. 

well as deep Ha imaging. The near-infrared photometric survey t^t^o i • i j- j • .i j j Tr> 

rlT — . o r . I j ,A7C- k J , , J ■ £ J /Tr.^ • • IC 348 was also extensively studied in the near- and mid-IR 

of| Lada&Lada| (|i995D, and later, deeper infrai-ed (IR) imaging Telescoue iLadaetall (l2006l) constructed 

studies were used to construct the luminosity function of the , ,. . ■rrnJ , ^ , , 

, ^ , J J ■ r . 11 J 1. . 11 optical/IR spectral energy distributions (SEDs) for all the cms 

cluster and provided information on the stellar and sub-stellar ^ i j- ■ a ^ ^ , i 



„ . ^ r; — r. , , „„o 177 1 — ^ i lUrvrvo r. -i • i tei" members and investigated both the frequency and nature of 

mass function (ILada et al.iil998t IMuench et aI.II2003c iPreibischi ^, ^ „ .• i w • i . irj . , . i 

bfifi'^ h the cir cumstellar disk population in the cluster. IMuench et al.l 

(l2007h performed a Spitzer census of IC348 and found that 



about 50 % of the known cl uster members possess circumstel- 

Send offprint requests to: B. Stelzer lar disks. ICurrie & KenvonI (12009.) presented new, deep MIPS 
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photometry of IC 348 that allowed a detailed characterization 
of the circumstellar disks around the known cluster members. 
Finally, millimetric observations directly revea led the presenc e 
and m a ss of 10 disks around IC 348 members (iLee et al.ll201 ih . 
iDahml (l2008h performed a spectroscopic investigation of accre- 
tion diagnostics for 40 near solar-mass members of IC 348 and 
derived accretion luminosities and rates for 14 of these stars. 

Based on all these studies, the basic properties of IC 348 can 
be summarized as follows: The mean age of the cluster mem- 
bers of 2 — 3 Myr is very interesting because it corresponds 
to the time where the structure of the disks of most young stellar 
objects changes from primordial, rather massive accretion disks 
to transitional and debris disks, and the point in time when plan- 
etary formation is thought to occur Furthermore, it offers in- 
teresting opportunities for comparisons to younger clusters like 
the Orion Nebula Cluster (1 ^2 Myr; Hillenbrand 1997) or the 
p Oph cluster (< IMyr; ILuhman & Riekelll999() . and older 
clusters and associ ations such as Upper Scorpius (~ 5 Myr; 
Ide Geus et al.]ll989h . 

The extinction of the cluster members ranges from Ay ^ 
1 mag up to Ay ^ 10 mag, with a mean value of ^ 3.5 mag. 
The majority of the cluster members is in the T Tauri stage of 
pre-main sequence evolution, and more than 70 stars can be clas- 
sified as "classical T Tauri stars" (CTTS) on the basis of their Ha 
emission. In the central parts of the cluster the active star forma- 
tion phase seems to be finished, and the cluster population repre- 
sents the outcome of a recent star formation event. About 10' to 
the south-west of the cluster center, however, a dense cloud core 
is found which contains several deeply embedded IR sources 
with extinctions exceeding ^ 20 mag in Ay as well as the 
very y oung molecular hydrogen jet HH 211 (i McCaughrean et aP 
11994 and the IC 348 MMS outflow (Eisl6ffeT etal.ll2003h . The 
star formation process is still underway in this region. 

IC348 is also well studied in the X-ray regime. Th e first 
X-ray observations of IC 348 by iPreibisch et al.l (Il996h were 
performed with the ROSAT observatory and led to the dis- 
covery of 116 X-ray sources. In September 2000, a deep X- 
ray image of IC348 was obtained with the Advanced CCD 
Imaging Spectrometer (ACIS) on board the Chandra X-Ray 
Observatory and led to the detection of 215 X-ray sources 
in the 17' x 17' field-of-view. These Ch andra observations 
were anal yzed by Preibisch & Zinneckej (1200 iL PZOl here- 
after) and IPreibisch & ZinneckeH (120021 PZ02 hereafter). In 
February 2003, IC 348 wa s observed with XMM-Newton 
( IPreibisch & Z innecke^ l2004l) . The XMM-Newton observation 
was strongly affected by solar particle flares and consequently 
the sensitivity was significantly reduced. Therefore, only a small 
number of new X-ray sources were detected in the field previ- 
ously observed with Chandra. Two more recent Chandra ob- 
servations ce ntered on the s o uth-w estern corner of IC 348 were 
analysed by iForbrich et alj (1201 ll) . They detected about one- 
third of the known cluster members in this area. 

While the analysis of these X-ray studies yielded numer- 
ous interesting results about the X-ray properties of the young 
stars in IC348, most of the above mentioned optical/IR stud- 
ies of the cluster were not yet available at the time. (The re- 
cent Chandra study of the fields in the south-west of the cluster 
by iForbrich et aP (iSol 1 ) focused on a comparison of radio and 
X-ray emission of YSOs without detailed consideration of the 
stellar properties.) The very substantial increase of optical/IR 
information about the cluster population during the last ^ five 
years clearly warrants a repeated effort to study the relations be- 
tween X-ray properties and stellar parameters. Another impor- 
tant aspect in this context is the availability of three new deep 



Chandra observations of IC 348, obtained in 2008 and spatially 
overlapping with the original field. The resulting Chandra 'mo- 
saic' comprises more than 90 % of the cluster population; see 
Sect. l4.1l for more details. 

During the last few years, several very deep and extensive 
X-ray studies of other low-mass star forming regions have been 
performed, which provide a good basis for comparisons of the 
X-ray properties of young stars at different ages and in differ- 
ent environments. The first of these studies was the Chandra 
Orion Ultradeep Project (COUP), a 10-day long observation 
of the Orion Nebul a Cluster (ONC) with Chandra/ ACIS (see 
iGetman et alj|2005 ), the deepest and longest X-ray observation 
ever made of a young stellar cluster A detailed investigation 
of the relation between the optical and X-ray properties of the 
~ 600 X-ray detected optically visible a nd well characte rized 
cluster members w as presented in Preibis ch et al.l (12005 ah and 
Stelzer et al. (2005). The XMM-Newton Extended Survey of the 
Taurus Molecular Cloud (XEST) covered the densest stellar pop- 
ulations in a 5 square degree region of the Taurus Molecular 
Cloud (see iGiidel et alJl2007 ') and provided X-ray data on 110 
optically well characterized young stars. 

The analysis presented in this paper uses all available 
Chandra observations of IC 348 to increase the spatial coverage 
of the cluster and the sensitivity in the region where the point- 
ings overlap. The X-ray source list is cross-coiTelated with an 
updated membership catalog that we have compiled from the 
literature. This yields stellar parameters, disk and accretion indi- 
cators and rotation rates for the majority of X-ray emitting YSOs 
in IC 348. We examine the dependence of X-ray emission level 
on these parameters and compare the results to X-ray studies of 
other nearby low-mass star forming regions such as Orion and 
Taurus with the scope to understand the origin of activity in pre- 
main sequence stars. 



2. The catalog of IC348 members 

2.1. The data base 

Since the studies by PZOl, the membership list of IC348 has 
changed significantly with more than 130 new members con- 
firmed by signatures of youth in optical low-resolution spec- 
tra and th eir position in the H R diagram above the 10 Myr 
isochrone ("Luhm an et aP I2003L henceforth LSM03). LSM03 
considered membership complete in the central 16' x 14' for 
M > 0.03 il/0 and Ay < 4 mag. However, additional mem- 
bers have been iden tified on the basis of IR spectroscopy 
(iLuhman et alj|2005bl henceforth LLM 05), and Spitzer mid-IR 
photometry ( Muench e t alJ l2007l) . of which about a dozen in 
the LSM03 completeness area. Several objects that had late-M 
spec tral types estimate d previously from HST near-IR photom- 
etry dNaiita et alJl2000l) have been rejected by LSM03 as cluster 
members on the basis of their optical spectra. 

Spitzer photometry from IRAC and MIPS for the members 
from LSM03 was discussed by Lada et al. (2006). They exam- 
ined the slope of the SED between 3.6. ..8 pm. for 263 objects 
with IRAC photometry in all bands, and distinguished proto- 
stars, disk-bearing stars and disk-less stars in four categories that 
we describe in Sect. 12.2! 

The present paper deals with the more than 90 % of known 
IC 348 members that are within the field of the presently avail- 
able Chandra observations. The catalog of optical/IR cluster 
members consists of the o bjects from Tab le 2 in LSM03, Table 2 
in LLM05 and Table 1 in lMuench et all ( 120071 ). These catalogs 
provide spectral types, effective temperatures (Toff), extinctions 
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Table 1 Definition of YSO class 



Class 


Q3-8 fiin 


I 


> -0.5 


II 


-O.5.... - 1.8 


II/III 


-1.8... -2.56 


III 


< -2.56 



Notes . Slope of the SEP fro m 3 — 8/im (adopted from lLada et alj 
ll2006h and lMuench etaP (l2007h ). 



(Aj or Ay), bolometric luminosities (£boi), and optical/near- 
IR photometry. For the objects from Luhman's papers Spitzer 
photometry is extracte d from Ladaetal. (2006) . A c cretio n 
rates are collected from lDahml (f2008i) : Muzerolle etaP (l2003h : 
[Mohantv et aL (20051) . and Ha equivalent widths are those listed 
in LSM03. 



2.2. Derived properties 

We have complemented the data base described above by several 
properties derived from the literature data. 

First, YSOs are classified in two different ways, (i) on the 
basis of the IR slope of their SED and (ii) on the basis of their 
Ha equivalent width. 

For the SED classifica tion the slopes a .s-sum derived by 
iLada et al.1 (|2006) and Mu ench et aP (l2007|) are used together 



with the boundaries defined by Lada et al 



2006*) that separate 

YSOs in different evolutionary phases (see Table [T]!. The four 
groups represent Class I objects (protostars). Class II (stars with 
thick disks). Class II/III (stars with "anemic" disks), and Class III 
(disk-less stars). The group termed "anemic" disks fills the gap 
between disk-bearing and disk-less stars in SpitzerflRAC color- 
color diagrams. This group shows a very small spread of [3.6]- 
[8.0] color which forms the lower envelope to thick disk stars. 

The distinction of classical and weak-line T Tauri stars 
(henceforth cTTS and wTTS) on the basis of the strength of Ha 
emission must consider the fact that the Ha equivalent width 
(Who) is not equivalent to the line flux as a result of its depen- 
dence on the continuum, which declines towards later spectral 
types. Therefore, when using Wna as an accretion indicator the 
threshold between cTTS and wTTS must be set at successively 
higher values for stars of later spectral types. We use thresholds 
of 3, 10, and 20 A for stars earlier than K5, between K6 and M3, 
and later than M3, respectively. This is similar to the criteria in- 
troduced by White & Ba sri (2003), who included an additional 
group with Who > 40A for stars later than M6. The omis- 
sion of this last step is justified by the separation of Class II and 
Class III sources in Fig. [T] 

A comparison of the two ways used to classify YSO evo- 
lution gives clues to the connection between the dissipation of 
dust and gas during the pre-main sequence. iLada et aTl (l2006t) 
found that in IC 348 there is a close relation between Ha equiv- 
alent width and IR excess. We confirm this in Fig. [T] using a 
more sophisticated division between accretors and non-accretors 
that depends on spectral type. The majority of Class II sources 
are classified as accretors on the basis of their Ha emission. 
Only one Class III source has an Ha equivalent width in the 
range typical for accretors, and only 6 out of 28 Class II/III 
sources. However, next to the majority of Class II/III sources, 
several Class II sources are classified as 'weak-lined' stars ac- 



cording to our criterion, i.e. they have disks but are not con- 
sidered to be actively accreting any more presumably because 
the gas has dispersed. These objects are likely the more evolved 
among the Class II and Class II/III stars. Another possible expla- 
nation for the presence of Class II sources below the cTTS cutoff 
line is variable accretion that would allow for some disk-bearing 
Class II and Class II/III objects to appear temporarily as wTTS. 
Stars with more than one Ha measurement (connected by ver- 
tical lines in Fig. [T]i show indeed, that Ha emission is strongly 
variable, and the equivalent widths measured at different epochs 
for some Class II objects do cross the boundary between cTTS 
and wTTS. For those stars we have used the highest value of 
Who as criterion for their cTTS/wTTS classification. 

It can also be seen from Fig.[T]that a few sources classified as 
Class I have not only a determined spectral type but also detected 
Ha emission, both unexpected for true protostars that ought to be 
invisible in the optical due to high extinction of their envelopes. 
In total, among the 15 Class I sources known in IC 348, a spectral 
type has been detected for ten of them and five have a measure- 
ment for Ha emission. The Ha equivalent widths of these stars 
are all above our 'accretor threshold'. These objects are prob- 
ably no genuine protostars but T Tauri stars with massive disks 
being respons i ble for thek shallow slope in the mid-IR. In fact, 
iMuench et al.l (l2007h state that the distinction between highly 
flared Class II disks and disk/envelope systems may require ob- 
servations at wavelengths longward of 10 fim. Such measure- 
ments are not available for the majority of our sample stars, and 
in any case they are not taken into account in the definition of 
YSO classes as described in Table [T] 

We have seeked to determine masses for all stars with tabu- 
lated Lboi and Tcff from an interpolation of evolutionary models. 
Different sets of models are known to be inconsistent betw een 
eac h oflier. LSM03 ha ve shown that the lBaraffe et all (1 19981) and 
Cha brier et al.l (l2000l) models yield the best description for the 
low-mass stars in IC348. These are also the only models that 
include very low mass stars (< 0.1 Mq) and extend into the BD 
regime. The HR diagram for IC 348 is shown in Fig. |2] where 
we anticipate the result from the cross-correlation with X-ray 
sources (see Sect. O. Our choice of models goes at the expense 
of an analysis of stars with Al > 1.2 Af0. However, as can be 
seen from Fig.|2]the vast majority of IC 348 members have lower 
mass. No mass could be determined for a few stars below the 
zero-age main-sequence. These are probably stars with edge-on 
disks seen in scattered light, consistent with their classification 
as YSOs with circumstellar matter (either Class II or Class I). 

3. X-ray data analysis 

IC348 was observed in September 2000 (ObsID 606, PI: 
Th. Preibisch) and three times in March 2008 (ObsIDs 8584 PI: 
N. Calvet, 8933 & 8944 PI: S. Wolk) with the Imaging Array 
of the Chandra Advanced CCD Imaging Spectrometer (ACIS- 
I), that provides a field of view of 17' x 17' on the sky. The 
aimpoint of observation 8584 was offset to the southwest with 
respect to pointing 606 and observations 8933 and 8944 were 
offset even more to the southwest; see Fig. |3] for a graphical 
representation of the Chandra observations. The total net expo- 
sure time of all four observations was 182 860 s (50.79 hours). 
Table |2]contains detailed information on the 4 observations. The 
roll angles (i.e. the orientation of the detector on the sky) are 
equal (288 - 289°) for the three observations from 2008. For ob- 
servation 606 the roll angle was 109.6° and hence the detector 
was rotated 180°. Therefore sources are only at the same detec- 
tor position for observations 8933 and 8944. Two of the CCDs of 
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log [K] 



Fig. 2 HR diagram for IC 348 compared to evolutionary tracks of iBaraffe et alJ (Il998h and IChabrier et alJ (l2000h . Isochrones are 
shown for 1, 2, 5, 10, 30, 100, 1000 Myr. The tracks are labeled in units of solar mass. X-ray detections are represented by filled 
circles and non-detections as crosses. IC348 members outside the Chandra fields are shown as open squares. The objects are 
color-coded according to their YSO class (see definition in Sect. l2.2l i. 



Table 2 Chandra observation log 



Obs.Id. 


Date Start [UT] 


Date End [UT] 


Exposure time 


Level 2 events 


Q!20oo [h:m:s] 


O2000 L ■ • J 


606 


2000-09-21 


19:58:42 


2009-09-22 11:09:29 


52 960 s 


347 817 


03:44:30.00 


-^32:08:00.0 


8584 


2008-03-15 


09:02:20 


2008-03-15 23:33:38 


50 140 s 


582 252 


03:44:13.20 


-^3 1:06:00.0 


8933 


2008-03-18 


17:35:07 


2008-03-19 15:12:52 


40 120 s 


294 144 


03:43:59.90 


-^31:58:21.7 


8944 


2008-03-13 


17:53:02 


2008-03-14 05:25:12 


38 640 s 


300 328 


03:43:59.90 


-^31:58:21.7 



the ACIS-S spectroscopic array were also turned on during the 
observations. However, since the PSF at the large off-axis angles 
at these detectors is strongly degraded, their point-source sensi- 
tivity is reduced; analysis was performed without using the data 
from ACIS-S. The basic data products of our observation are the 
four Level 2 processed event lists provided by the pipeline pro- 
cessing at the Chandra X-ray Center, that list the arrival time, 
location on the detector and energy for each of the 1524541 
detected X-ray photons. We combined the four pointings with 
the mergejill script, a Chandra contributed software that makes 
use of standard tools. The merged events file is shown in 

Fig. [3] Overlaid are crosses at the positions of the IC 348 mem- 
bers showing the spatial distribution of the YSOs. 

3. 1. Source detection and X-ray source catalog 

The source detection was caiTied out in a two-step process. The 
first detection step was performed in a rather aggressive man- 

' Chandra Interactive Analysis of Observations, version 4.2: 
|http://cxc .harvard.edu/ciao/index.html| 



ner in order to find even the weakest possible sources, deliber- 
ately accepting some degree of false detections. In the second 
step, this list of potential sources was then cleaned from spuri- 
ous detections by a detailed individual analysis. We employed 
the WAVDETECT algorithm (a CMO mexican-hat wavelet source 
detection tool; Freeman et all l2002h for locating X-ray sources 
in our merged image, and used a detection threshold of 10^^ and 
wavelet scales between 1 and 16 pixels. The result was a list of 
372 prospective sources. Three sources were added, which have 
been detected with WAVDETECT in an XMM-Newton observa- 
tion of IC 348. Another 17 possible sources were added, which 
were identified by cross-correlating UKIRT observations with 
the Chandra image by eye. These 20 sources were first added 
as potential sources to the sourcelist for ACIS Extract. ACIS 
Extract then constructs an extraction region at the position of 
these sources and estimates the probability for being a source. 
The result was a final catalog with 393 potential X-ray sources. 
To clean this catalog from spurious sources, we performed a de- 
tailed analysis of each individual candidate source with the ACIS 
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3.4 



3.5 



3 6 

logLJK] 



3.7 



3.8 



Fig. 1 Equivalent width of Ha vs. effective temperature. The 
color-codes refer to the Spitzer YSO classes: Class III (green; 
light grey in b/w print), Class II/III (blue; dark grey in b/w print), 
Class II (red; medium grey in b/w print) and Class I (yellow 
and larger). Different mass ranges are symbolized with different 
plotting symbols: 0.1 - 0.25 Mq (star symbols), 0.25 - 0.6 Mq 
(circles), 0.6 — 1.2 M0 (triangles) and all others (squares). 
Vertical lines connect different epochs of Wua measurements 
for individual stars. The black line separates accreting cTTS 
(top) from non-accreting wTTS (bottom). The cross on the black 
line denot es spectral type M6. Fo r objects with Toff smaller than 
this value IWhite & Basril (l2003h have placed the dividing line 
between cTTS and wTTS at 40 A. However, our boundary of 
20 A is in better agreement with the separation of Class II and 
Class III sources. 



Extract (AE hereafter) software packag^ dBroos et alj|2010h . 
A full description of t h e procedures used in A E c an be found 



in 



Getman et al.1 (l2005l) . rTownslev et alj (l2003h . and lBroos et"al] 
(I2007h . The following three steps were performed by AE in or- 
der to prune our input catalog from spurious detections (includ- 
ing afterglows): 

1 . Extraction regions were defined such as to include 90% of 
the source photons in the PSF (or less in the case of other 
nearby sources), and source events were extracted from this 
region. ARFs and RMFs are constructed for each source, and 
the ARFs are corrected for the light missed by the finite ex- 
traction regions. 

2. Local background events were extracted after masking all 
the sources in the catalog. 

3. The Poisson probability {Pb) associated with the "null hy- 
pothesis", i.e. that no source exists and the extracted events 
are solely due to Poisson fluctuations in the local back- 
ground, was computed for each source. 

After 6 iterations of this pruning procedure our final cata- 
log consisted of 290 X-ray sources, for which the probability 



of being a background fluctuation is less than 0.1 %. For 287 
sources of those the probability is less than 0.01 % and 266 
sorces have a very low probability of 10^^. Some of the ba- 
sic source properties derived by the AE software, such as the 
net (i.e., background-subtracted) counts in various energy bands, 
the median photon energy (ii^med), statistical test for variabil- 
ity, and a measure of the incident photon flux, are reported in 
Table I?] Sources are sorted by increasing right ascension and 
identified by their sequence number (Col. 1) or their lAU desig- 
nation (Col. 2). 

3.2. Calculation of X-ray fluxes 

An accurate determination of the intrinsic X-ray source lu- 
minosities requires good knowledge of the X-ray spectrum. 
However, for the majority of the X-ray sources the number of 
detected photons is too low (less than 100 net counts) for a de- 
tailed spectral analysis. 

An estimate of the intrinsic, i.e. extinction corrected. X-ray 
luminosity for sources that are too weak for a detailed spectral 
an alysis can be obtained with the XPHOT softwarfl developed 
bv lGetman et alJ (l2010l) . XPHOT is based on a non-parametric 
method for the calculation of fluxes and absorbing X-ray col- 
umn densities of weak X-ray sources. Our procedure is based 
on a grid of simulated X-ray spectra with the following param- 
eters: two temperatures, the emission measure ratio of the two 
spectral components, and the absorbing column. The X-ray ex- 
tinction and intrinsic flux of a given X-ray source are estimated 
from a sophisticated comparison scheme between the observed 
median energy and the observed source flux and analogous val- 
ues obtained from the simulated spectra. This method requires at 
least 4 net counts per source (in order to determine a meaningful 
value for the median energy) and can thus be applied to 250 of 
our 290 sources. 

An estimate of the observed (i.e. not the intrinsic) X-ray 
flux is also computed by AE. This quantity, called FLUX2 and 
given in units of photons/s/cm^, is calculated from the number 
of detected photons and using a mean value of the instrumental 
effective area (through the Ancillary Response Function, ARE) 
over energy. It should be noted that the FLUX2 values suffer 
from a systematic uncertainty with respect to the true incident 
flux, because the use of a mean ARE is only correct in the hy- 
pothetical case of a flat incident spectrum, an assumption that is 
probably not fulfilled. Nevertheless, the FLUX2 represents the 
best estimate for the photon flux that can be obtained for com- 
paratively weak sources. The product of FLUX2 and i?nicd pro- 
vides an estimate for the energy flux of the source. When com- 
pared to the fluxes derived with XPHOT a rather good linear 
correlation is found but the product FLUX2 ■ Ej^cd results in 
systematically higher values. We calibrated the FLUX2 values 
using the correlation with the XPHOT fluxes to obtain source 
fluxes for the 40 sources that are too faint for XPHOT. 

This way, the X-ray flux of all 290 detected sources has 
been determined in a consistent way. We cross-correlate the X- 
ray source list with the catalog of IC 348 members from Sect. |2] 
using a match radius of 1.5". Larger cross-correlation radii did 
not substantially increase the number of matches, e.g. a radius 
of 3.0" yields only four additional identifications, an increase 
of the X-ray sample for IC 348 members of about 2 %. We are, 
therefore, confident that nearly all X-ray emitting IC 348 mem- 
bers have been recognized by our procedure. The flux of all X- 
ray detected IC 348 members is given in Table [3j A distance of 



http://www.astro.psu.edu/xray/docs/TARA/ae_users_guide.html 



http://www.astro.psu.edu/users/gkosta/XPHOT/ 
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Table 3 X-ray fluxes and upper limits 



I h ft M ri vff \ Fl 

y^llLlllLlf Ll LI-/ 




lop- f flap'' 
[erg/cm /s] 




1 
i 


^ 1/1 Q/1 
< — 14.o4 


Ci'lAA'i.'^ Q<_i_T7inn/1 K 
Uj44j J. J j+jZ1UU4.Z> 


z 


— iz.oz 






— iZ.iU 




/I 


— lo. ( 1 




c 

J 


— iZ.o4 


U j44 JO. V J+ jZUd4 J .4 


D 


1 o no 
— iz.uy 




7 


< —14.35 




8 


< -14.07 


034439.15+320918.1 


9 


-12.25 


034424.66+321014.8 


10 


-14.28 


034507.96+320401.8 


11 


-12.32 




12 


< -14.78 




12 


< -14.63 


034359.67+320154.1 


13 


-12.53 



Notes. '"^ Identifier from LSM03, LLM05 or lMuench et"an d2007h for 
all IC 348 members observed with Chandra. '''' Flag is set if the flux 
was obtained from AE (FLUX2). 



SIOdc (lHerbigll998h is assumed for conversion of the flux to lu- 
minosity. The resulting intrinsic X-ray luminosities range from 
1028.38 to 1031-93 erg s-1. 



Table 4 X-ray detection statistics for different YSO classes in 
IC348 



Mass Range [JI/q] 




I 


n 


Will 


III 


All 


all IC 348 members 


iVul 


8 
6 
0.57 


44 
61 

0.42 


36 
32 
0.53 


94 
27 
0.78 


187 
129 
0.59 


M < 0.1 


iVdct 
iVul 



1 


1 

19 
0.05 


1 

14 
0.07 


3 
7 
0.30 


5 

41 

0.11 


0.1 < M < 0.25 


A^dct 

N 



2 


10 
28 
0.26 


17 
15 
0.53 


23 
13 
0.64 


50 
58 
0.46 


0.25 <M< 0.6 


A^dct 
iVul 

Ndat 

N 






11 
9 
0.55 


9 
1 

0.90 


30 
4 
0.88 


51 
15 
0.77 


0.6 < M < 1.2 


iVul 
JV 


4 
1 

0.80 


14 
2 
0.88 


6 

1.00 


21 

1.00 


45 

3 

0.94 



3.3. Determination of upper limits 

For th ose 129 stars from LSM03, LLM05 and ' Muench et al] 
(|2()07!), which had no counterpart in our X-ray source catalog, 
upper limits had to be determined. We did this using again the 
AE software package. Extraction regions were defined the same 
way as the 90% contours of the local PSF to determine the num- 
ber of counts in the target aperture and an estimate of the local 
background. To calculate upper limits for t he received sourc e 
counts from AE we used the program bayes.f dKraft et alJl99lh . 
This program performs the calculation of the upper limit on a 
Poisson process with background, incorporating the effects of 
uncertainty on the background and signal acceptance. Upper 
limits for the number of source counts were determined for a 
confidence level of 0.9. For each object we divided the upper 
limit for the source counts by the local exposure time. To ob- 
tain upper limits for the luminosity we determined a conver- 
sion factor (CF) between count rate and luminosity. As coun- 
trate and XPHOT luminosity are correlated for T Tauri stars, 
this was done with a sample of 180 stars with known spec- 
tral type from our AE source catalog. The mean of the CF is 
3.74462 • lO^-^ erg/counts. The CF depends essentially on ab- 
sorption, and we verified that the sample of detected stars from 
which it was calculated has a similar mean Aj value as the sam- 
ple of undetected stars to which the CF was applied. The upper 
limit fluxes for those latter ones are listed in Table [3] together 
with the fluxes of the detected IC 348 members. 



4. Results 

4. 1. X-ray detection statistics 

Our combined membership list comprises 345 YSOs from spec- 
tral type B5 to M9. Only 29 of them are located outside the 
sky region observed with Chandra. We have X-ray detected 187 
YSOs (59 % of the observed sample); see Table |4] Many of the 
new members identified after the study of PZOl are too faint for 



detection in the existing Chandra data, even after combining the 
original image (Obs-ID 606) with the more recent exposures. 

It is evident from the position of the 129 X-ray undetected 
IC 348 members in the HR diagram (Fig.|2]i that these are mostly 
very-low mass stars and BDs. A closer look at the X-ray prop- 
erties of the lowest mass objects in IC 348 is shown in Fig. |4] 
(left) where we display X-ray luminosities for spectral class M. 
Only 4 of the 39 known presumable BDs (defined here as objects 
with spectral type M6 and later) within the Chandra FOV are de- 
tected. The coolest X-ray detected IC 348 member has spectral 
type M7.5. PZOl claimed the X-ray detection of 7 BDs and BD 
candidates. However, only two of them were detected by their 
automatic source detection routine and the other five by collect- 
ing the counts within a 3" region around the optical position. 
In our new analysis we have used a more sophisticated source 
detection procedure. Moreover, the updates of the stellar param- 
eters in the recent IC 348 membership studies have shifted two 
of the objects previously considered BDs to earher spectral type, 
i.e. stellar mass (see LSM03). 

4.2. X-rays and evolutionary stage 

Fig. |5] shows color-color diagrams for Spitzer/IRAC and 
SpitzerMTPS data for Lada Class I, II, Will, and III sources. 
X-ray detected and undetected stars are represented by filled 
and open symbols respectively. In the left panel we overplot 
the Class II areas defi ned by Megeath et al. (20 04) and by 
iHartman n et al.l (l2005l) . As was shown already by iLada et aD 
(2006), the YSO classification based on the SED slope corre- 
sponds to distinct regions in the Spitzer color-color diagrams, 
in particular those involving [3.6] - [8.0]. This classification is 
in good agreement with the color-color selection for Class II 
sources introduced by Megeath et al. (2004), although the latter 
one includes a significant number of stars in t he Class II area tha t 
are considered transition disks acco rding to Lada et al. ( 20061) . 
The Class II criteria introduced bv IHartmann et aL ( 20051) are 
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more conservative and include almost exclusively objects with 
Lada Class II flag but they miss a good fraction of them. 

The X-ray detection statistics for the different YSO classes 
in IC 348 are summarized in Table |4] Columns 'N^et and '-/Vuf 
are the number of stars detected and undetected ('upper limits') 
with Chandra, respectively. In the last column we give the to- 
tal number of stars, regardless of YSO class. Only 8 confirmed 
IC 348 members within the area surveyed by Chandra have no 
YSO classification. The detection fraction is highest for the most 
evolved stars. Class III, regardless of mass. For each YSO class 
there is a clear trend of increasing detection fraction with in- 
creasing mass. In particular, for stars with M > 0.6 Mq all 
Class II/III and III sources are detected, and the total detection 
statistic regardless of YSO class is > 94 %. On the other hand 
only ^ 10% of very low-mass objects with M < 0.1 Mq are 
X-ray detected. 

LSM03 consider membership complete in the central 16' x 
14' for M > 0.03 Mq and Ay < 4 mag, but this estimate might 
be biased at the low-mass end against disk-less or strongly ab- 
sorbed members as corroborated by the identification of further 
cluster members with IR spectroscopy and Spitzer mid-IR pho- 
tometry, several of which in t he completeness area of LSM 03 
(LLM05: lMuench et al] (12007')). The Spitzer color selection ex- 
cludes by definition the identification of Class III sources. 

The X-ray luminosities of the different YSO classes are com- 
pared in Sect. 14.41 Here we point at -Emcd as a rough measure for 
spectral hardness and/or extinction. The median energy is dis- 
played for X-ray detected M stars of IC 348 in the right panel 
of Fig. m For clarity only the spectral class M is displayed. The 
mean of the median energy ((i^med)) for the full sample irre- 
spective of spectral type is shown for each YSO class at the right 
end of the graph together with its standard deviation. It strikes 
that (i?mcd) of Class I objects is much higher than that of all 
other YSO classes. This is likely due to high extinction from 
protostellar envelopes or thick disks absorbing the soft photons. 




-2-10 1 

log (UL,j 

Fig. 6 X-ray versus bolometric luminosity (a small number 
of objects, not relevant for this study, falls outside the plot- 
ted Lboi range); dotted lines represent constant Lx/iboi ra- 
tios 10^'^, 10^'* and 10^^; the dashed line is the linear regres- 
sion fit derived by Preibischet al. (2005a) for the ONC corre- 
sponding to log(ix/iboi) ~ —3.7; the solid lines are linear 
regression fits to the IC348 members in the luminosity range 
log (L/Lq) — —1.5... + 0. 5, which correspond s roughly to the 
upper mass boundary of the lBaraffe et al.l (Il998l) models: black, 
red and green are for the whole sample, only Class II and only 
Class III sources, respectively. 



4.3. X-ray and bolometric luminosity 



The ratio between X-ray and bolometric luminosity is thought 
to be a crucial indicator for the level of magnetic activity in 
young stars. From X-ray surveys in various star forming re- 
gions it has emerged that the characteristic emission level of 
most low-mass pre-main sequence stars is concentrated between 
log(Lx/Lboi) ~ —3 and —4 but with a spread of at least 
3 dex, and there are indications for a sligh t decline of the 



log (Lx/^boi) level with de cr easing mass (e.g. Flaccomio et alj 



l2003tlPreibisch et al.ll2005allbl:lGrosso et alJl2007l) 



4.3.1 . Sample without distinction of YSO class 

Fig.|6]shows the vs. Lboi diagram for IC 348. The sample of 
known IC 348 members extends to lower bolometric luminosi- 
ties and masses than that of Orion, possibly due to both environ- 
mental and observational effects (e.g. higher extinction in Orion) 
and the closer distance. We carried out a linear regression fit for 
the range lo g,(L/ Lq) = —1. 5... + 0.5 using the EM algorithm 
in ASURV (iFeTgelson & Nelsonill985h . Our choice for the lu- 
minosity cutoffs are motivated by the properties of the IC 348 
sample (avoiding at the high end X- ray faint mass i ve sta rs out- 
side the mass range covered by the ' Baraffe et alJ (Il998l) mod- 
els and avoiding to be dominated by upper limits in X-rays at 
the low end). The resu lt is overplotted in Fig. |6] together with 
the bestfit obtained by iPreibisch et all (12005 a) for the ONC in 



a similar range of bolometric luminosities. The ONC study was 
performed for stars with log (L/Lq) < 1 and without lower cut- 
off to the bolometric luminosity. However, effectively the ONC 
sample extends down to log(L/io) — —1.5 , i.e . the same 
as our low-luminosity cutoff. With respect to IPreibisch et al] 
(i2005ah . we have chosen a lower value for the luminosity cut- 
off at the high end because at the age of IC 348 stars with 
log {L/Lq) > 0.5 are mostly radiative (cf. Fig. |2]l and, there- 
fore, their X-ray emission - if any - may not have the same ori- 
gin as in lower-mass stars. 

The slope we determined for IC 348 in the luminosity range 
described above (1.27 ± 0.08) is significantly steeper than the 
one measured from COUP for t he ONC (1.04 ± 0.06 ) and from 
XEST for Taurus (1.05 ± 0.06; iTelleschi et ani2007 ). While in 
Orion and Taurus the linear fit of L^ vs. Lboi is consistent with 
a constant Lx/Lboi ratio, in IC348 we observe a decrease of 
ix/iboi with decreasing bolometric luminosity. A possible rea- 
son may be the particular distribution of YSO types and their 
different X-ray luminosities. Note, however, that the COUP and 
XEST samples we refer to are not directly equivalent to our sam- 
ples in IC348 because of diffe rent ways the subsamples were 
defined. IPreibisch et al.l (l2005ah distinguish accretors from non- 
accre tors on the basi s of 8542 A calcium emission. Similarly, 
Tell eschi et al.l (l2007l) separate accretors from non-accretors on 
the basis of Ha emission, while we have chosen to work with 
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the Spitzer YSO classification because it yields the largest sam- 
ple size. 

4.3.2. Class II vs. Class III and comparison with the ONC 

When the same analysis is carried out separately for the subsam- 
ples of Class II and Class III sources in IC 348 we find slopes 
of 1.65 ± 0.22 and 1.13 ± 0.11, respectively. These results are 
graphically demonstrated in Fig.|6l To summarize, we find that 

- the Lx — ^boi relation of Class III objects in IC 348 is in 
agreement with that of the non-accretors in the ONC and is 
roughly constant across the examined range of bolometric 
luminosities; 

- the Class II sources in IC 348 show a steeper slope, i.e. de- 
creasing Lx — iboi ratio towards fainter stars. IC 348 hosts 
a large number of Class II sources with low luminosities (i.e. 
low masses) - see e.g. Table |4]- and these objects have low 
X-ray luminosities or they are not detected. 

Only recently YSO classifications for ONC stars based on 
SpitzerflRAC photometry have become available. We extracted 
these data from the YSOVAR databas^ In Fig. Qwe compare 
the X-ray vs. bolometric luminosities for IC 348 with those of 
the ONC separately for Class II and Class III sources. 

According to Fig. iTjleft), in contrast to IC 348 there is no 
evidence for an X-ray deficiency of low-luminosity diskbearing 
stars in the ONC. However, a look into the HR diagram for the 
ONC, presented in Fig.[8j suggests that the bolometric luminosi- 
ties of these stars (i.e. the ones with logLboi ^ —1... — 1.5) 
may be underestimated. This group is characterized by an ob- 
vious excess of Class II sources with isochronal age around 
lOMyr or more in the HR diagram, in contradiction with the 
general paradigm that Class II stars are on average less evolved 
than Class III stars as well as with plausible timescales for disk 
dissipation of just a few Myr. We hypothesize that many of 
these low-luminosity, low-mass Class II objects are either ob- 
served in scattered light, the central star being fully obscured by 
an edge-on disk, or partially obscured by disk warps or other 
circumstellar structures the so-called AATau phenomenon. A 
number of AATau-like stars in the ONC were identified by 
iMorales-Calderon et a n (l2011h and are marked in Fig.[8]by large 
annuli. Not all of the cluster members with untypically low 
bolometric luminosity for their effective temperature display the 
AA Tau features but this is not in contradiction with our hypothe- 
sis because occultations by disk warps are transient phenomena. 
In fact , these objects were identified bv 'Morale s-Calderon et alj 
(1201 Ih on the basis of their time-variability. Moreover, some of 
these objects may be permanently obscured systems with edge- 
on disk. Both of the described scenarios would imply that the 
extinction correction for these stars has been significantly un- 
derestimated, thus placing the true position of the stars in the 

— iboi diagram further to the right and in better agreement 
with the trend that we see in IC 348. The fact that no such Class II 
stars with low Lboi but high are present in IC 348 might be 
explained by the older age with respect to the ONC, going along 
with the dispersal of the circumstellar material that provides the 
extinction. 

We pass now to the consideration of the C lass III sources 
(Fig. Illright). In contrast to the analysis of jPreibisch et alJ 
(2005a) for non-accretors, the ONC Class III sample shows evi- 
dence for a steeper — Lboi relation than the one correspond- 
ing to a constant luminosity ratio. This is in line with the result 

http://ysovar.ipac.caltech.edu 




3.7 3.5 3.5 

Log T.„ [K] 



Fig. 8 HR diagra m for the optical COUP sample on the 
ISiess et aL 1 (l2000h models. Class III ( green; light grey in b/w 
print) and Class II (red; dark grey in b/w print), AA Tau like stars 
(large blue annuli) and IC 348 members without YSO classifica- 
tion (small black points). 



for Class II sources considering the hypothesis made above con- 
cerning AATau like stars. We could not detect this deficiency 
of low-luminosity X-ray faint Class III sources in IC 348 but this 
may be due to the limited sample size. From Fig.|7]it appears that 
the IC 348 and ONC Class III distributions are rather similar. 

4.4. X-ray luminosity functions 

Cumulative distributions of the X-ray luminosities are tradition- 
ally used to compare subsamples of pre-main sequence stars. In 
the past, the X-ray luminosity functions (XLF) derived for var- 
ious star forming regions have been vigourously debated, as it 
could not be established if circumstellar matter influences these 
distributions. Most studies that have distinguished accreting and 
non-accreting stars on basis of the strength of their Ha emission 
found that c TTS have systematically lower X-ray luminositie s 
than wTTS (IStelzer & NeuhauseillloOll: iFlaccomio et al.ll2003h . 
On the other hand, no differences were found between stars with 
and without disks as diagnosed by IR excess (e.g. PZOl). 

We have constructed cumulative distributions of both the X- 
ray luminosity (i.e. the XLFs) and the Lx/-^boi ratio. We com- 
puted XLFs for YSOs of different SED classes (to examine the 
influence of disks) and for different strength of Ha emission 
(to examine the influence of accretion). To avoid biases related 
to the known dependence of on stellar mass, the analysis 
was carried out in three mass bins chosen such that they sam- 
ple adequately the low-mass IC 348 population: 0.1 — 0.25 M©, 
0.25 — 0.6 A/0, 0.6 — 1.2 Mq. The significance of the results is 
verified with two-sample tests that yield the probability for re- 
jecting the hypothesis that the observed luminosity distributions 
are drawn from the same parent distributions. Both the cumula- 
tive distributions and the two-sample tests are computed within 
the ASURV environment. The calculations take into account the 
upper limits for undetected stars. 

The results are displayed in Figs. |9] and [TO] In the left hand 
panels we compare the XLFs of Class II, llflll, and III sources. 
Class I objects are not considered because of poor statistics. 
There is a pronounced difference in the X-ray luminosities of 
the three YSO classes for the lowest mass bin (0.1 — 0.25 Mq), 
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Table 5 Statistics for detections, median of the logarithm of X-ray luminosity and Lx/^boi ratio, and 2-sample tests of different 
YSO classes 

















2S-Test 


2S-Test 


2S-Test 




Class III 


Class II/III 




Class II 


Class III -II 


Class III -II/III 


Class II/III -II 


M[Me] 


N/Nui 


log Lx 


N/Nui 


log Lx 


iV/iVui 


log Lx 


Prob.[%] 


Prob.[%] 


Prob.[%] 


0.1 - 0.25 


36/13 


28.76 


32/15 


28.51 


38/28 


28.47 


> 99 


> 97 


< 90 


0.25-0.6 


34/4 


29.20 


10/1 


29.43 


20/9 


29.04 


< 90" 


< 90 


< 90" 


0.6-1.2 


21/0 


29.55 


6/0 


29.90 


16/0 


29.76 


< 90 


< 90 


< 90 


M[MqJ 


iV/A^ui 


logLx/Lboi 


iV/iVui 


log Lx/Lboi 




log Lx/Lboi 


Prob.[%] 


Prob.[%] 


Prob.[%] 


0.1-0.25 


36/13 


-3.78 


32/15 


-3.87 


38/28 


-4.18 


> 99 


< 90 


> 97 


0.25-0.6 


34/4 


-3.60 


10/1 


-3.74 


20/9 


-4.15 


< 90'' 


< 90 


< 90 


0.6-1.2 


21/0 


-3.59 


6/0 


-3.27 


16/0 


-3.58 


< 90 


< 90 


< 90 



Notes. The logrank test gives probabilities of > 97 % and > 93 % for Class III vs. Class II and Class II/III vs. Class II. The logrank test 
gives a probability of > 99 % for Class II vs. Class II. 



where Class III sources have higher Lx than both Class II/III and 
Class II and the latter two classes are indistinguishable. This is 
also the mass range with the highest number statistics compris- 
ing 106 stars. However, these low-mass samples have a large 
number of upper limits. A marginal difference between Class II 
objects and the other two YSO groups is found for the interme- 
diate mass bin (0.25 — 0.6 Mq) only in the logrank test while 
all other types of two-sample tests could not detect a differ- 
ence with better than 90 % probability. For the highest mass bin 
(0.6 - 1.2 Mq) the XLFs of the three types of YSOs can not be 
distinguished. The graphical results are supported by the two- 
sample tests. The statistics are summarized in Table|5]where 'N' 
is the total number of stars and 'N^^^ is the number of undetected 
stars. The luminosities refer to the median of each sample, and 
the numbers in the last three columns represent probabilities for 
two distributions being different obtained from two-sample tests. 

Also in terms of Lx/iboi there is a marked difference be- 
tween the three YSO classes only in the lowest mass bin. In con- 
trast to the XLFs, where Class II/III objects behave similar to 
the Class II sources, in Lx/Lboi the Class II/III objects are more 
similar to the Class III sources (compare uppermost left panels 
in Figs. |9] and [Toll. Again only the logrank test gives Class II ob- 
jects different from Class III for the intermediate mass bin. All 
other distributions of Lx/Lboi for higher masses are undistin- 
guishable. These results confirm our findings from Fig. |6] that 
the Lx — iboi relation in IC 348 is steeper for Class II than for 
Class III sources. 

An analogous investigation was performed for cTTS and 
wTTS, discriminated by the equivalent width of Ha as described 
in Sect. 12.21 The XLFs and the cumulative distributions for 
L^/ Lboi for cTTS and wTTS are shown in the right panels of 
Figs.|9]and[T0]for the same mass bins that we have examined for 
the samples classified by their IR excess. Again, a significant dif- 
ference between the groups is found for the lowest mass bin with 
wTTS being more X-ray luminous than cTTS; see Table |6] for a 
summary of the results. In the intermediate mass bin, where the 
Class II stars were different from the other two samples only in 
the logrank test, the cTTS and the wTTS are different at high sig- 
nificance for all two-sample tests if is considered but, again, 
only in the logrank test if L^/ Lboi is considered. 

Note, that the IR excess and the Ha emission samples are 
not identical and their comparison is not straightforward. The 
number of stars with measured aa-g^m is larger than those 
with Ha data. For the lowest mass bin, most of the Class II/III 
sources belong to the wTTS group. However, there are also sev- 
eral Class II sources that are classified as non-accretors. We have 



Table 6 Statistics for detections, median of the logarithm of X- 
ray luminosity and Lx/Lboi ratio, and two-sample tests for cTTS 
and wTTS 



wTTS 



cTTS 



2S-Test 
wTTS - cTTS 



M[Mq] 


N/N^i 


log Lx 


N/N^i 


log Lx 


Prob.[%] 


0.1-0.25 


55/18 


28.74 


26/20 


28.18 


> 99 


0.25 - 0.6 


38/3 


29.38 


8/4 


28.51 


> 95 


0.6 - 1.2 


16/0 


29.76 


15/0 


29.78 


< 90 


M [Mq] 


N/Nui 


log 7^ 


N/Nui 




Prob.[%] 


0.1-0.25 


55/18 


-3.78 


26/20 


-4.09 


> 99 


0.25 - 0.6 


38/3 


-3.59 


8/4 


-4.57 


< 90" 


0.6 - 1.2 


16/0 


-3.32 


15/0 


-3.64 


< 90 



Notes. The logrank test gives a probability of > 96 %. 



shown that the trend of a mass dependence of the X-ray emission 
level throughout the various evolutionary stages of YSOs is ro- 
bust against details of the sample definition. 

4.5. X-ray luminosity and mass accretion 

In the past, differences between the X-ray luminosities of cTTS 
and wTTS have been ascribed to a possible connection between 
X-ray emission and accretion. The interpretations have involved 
opposite views: In one scenario X-ray emission is suppressed by 
accretion, e.g. by a direct influence of accretion onto the dynamo 
efficiency through reduced convection or through reduced coro- 
nal heating as a result of the higher plasma d ensities involved 
in the accretion columns (e.g. iPreibisch et al. 2005a) or by oc- 
culting the X-ray emitting corona (lOregorv et al.ll2007h . In the 
alternative picture X-ray emission determines whether accretion 
takes place or not through photoevaporation of disks. This lat- 
ter hypothesis has been bolstered by an empirical anticorrelation 
of X-ray luminosity and mass accretion rate observed in COUP 
(iDrake et al. 2009). 

Only 18 IC 348 members have measured mass accretion rate. 
In Fig. [TT] we present the relation between Lx and Mace for 
IC 348. Note, that at least two stars in IC 348 with 'anemic' disks 
are accreting. The top panel of this figure shows that the appar- 
ent positive correlation between and ilfacc is a result of the 
well-kn own dependence o f mass accretion rate on stellar mass 
(see e.g. lNatta et al.ll2006h . Three of the objects in the group of 
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Fig. 1 1 X-ray luminosity and mass versus mass accretion rate. 
Colors define YSO classes as in Fig.[T] 



strong accretors (empirical dividing line of log il/acc = —9) have 
no mass as signed becaus e they are outside the mass range cov- 
ered by the iBaraffe etal] d 1 998i) models (Af > 1.2 Mq). All in 
all, however, all stars in this group have similar mass but span 
a range of two decades in accretion rates and more than one 
decade in X-ray luminosity. The Spearman's test does not yield 
evidence for a correlation between the two parameters for the 
IC 348 members with strong accretion. However, the small sam- 
ple size precludes any firm conclusions. 



5. Discussion and conclusions 

Our combined analysis of all existing Chandra observations pro- 
vides the so far deepest X-ray survey for IC 348. It is difficult to 
give a number for the sensitivity limit due to the partial overlap 
of the available images. The exposure time for the deepest part 
of the merged mosaic is 180 ksec and the faintest X-ray detected 
cluster member has a luminosity of logix [crg/s] = 28.1 Our 
detection rate for BDs in IC 348 is 10 %, in an area surveyed by 
Chandra that misses only 2 out of 41 substellar objects. 

The magnetic activity at and beyond the low-mass end of the 
main-sequence is generally poorly constrained due to a lack of 
sensitive enough observations in the X-ray regime. Even COUP, 
the so far deepest X-ray image of a star formation region, has 
yielded detections for only 26 % of the known BDs, presum- 
ably due to the considerable extinction in the ONC. We caution 
that the COUP BD sample may be contaminated by a few ob- 
jects with mass above the substellar limit due to the fact that the 
BD classification was based mostly on near-IR spectra that are 
known to yield systematically later spectral types compared to 
optical spectra (see Preibisch et al. 2005b, for more details on 
the sample definition). Nine B Ds were detected in Taurus during 
the XEST dGrosso et alJ2007l) using the same definition for BDs 



as we do in this paper for IC 348 (spectral type M6 and later). 
The XEST sample comprised 17 BDs known at the time, i.e. the 
X-ray detection rate was about 50 %. It is important to note that 
many of the X-ray detected BDs in Taurus have untypically high 
bolo metric luminos i ties pl acing them above the 1 Myr isochrone 
of thelBaraffe et akld 19981) models (see e.g. Fig. 2 in lGrosso et al.l 
'2007). As a consequence, for a given Lx/^boi level their X-ray 
luminosity is higher than for IC 348 BDs explaining the higher 
detection rate. 

We have examined cumulative distributions of absolute and 
fractional X-ray luminosity in three mass ranges: 0.1 — 0.25, 
0.25 — 0.6, and 0.6 — 1.2 A/q. The choice of these mass bins was 
driven by the requirement to have sufficient stars in each sub- 
sample for a statistical analysis. The well-known increase of the 
median X-ray luminosity with increasing mass underlines once 
more that the discrimination of different mass bins is manda- 
tory. However, the details of the derived XLFs also depend on 
the choice of the boundaries for the mass bins, and more gener- 
ally on the pre-main sequence models. For this reason, results for 
different star forming regions are difficult to compare. In COUP, 
e.g., the Siess et al. (2000) models were used, and the mass bins 
that were chosen are different from ours motivated by the high 
number of stars with solar mass and above. 

In addition to the issue with stellar masses, the literature 
provides differing definitions for characterizing the evolution- 
ary stage of pre-main sequence stars. In COUP, e.g., disks 
were characterized by near-IR excess effectively biased towards 
warmer disks with respect to our Spitzer YSO classification. 
Other COUP subsamples represent accreting and non-accreting 
stars separated by the strength of calcium emission. In XEST 
YSOs were distinguished into accretors and non-accretors on 
the basis of Ha emission while no disk-status was assigned to 
the stars. We have examined for IC 348 both the influence of the 
presence of disks and of ongoing accretion on the X-ray emis- 
sion level. Our results hold for both types of classifications. In 
particular, we find that differences exist between groups of YSOs 
in different evolutionary stages at low masses while they disap- 
pear at higher mass. This confirms a trend that was seen in COUP 
(Preibisch et al. 2005a). An opposite tendency was claimed in 
XEST: different XLFs for cTTS and wTTS at M > 0.3 Mq but 
not for lower masses ( Telle schi et al ] |2007.) . This result does not 
hold for IC 348 according to our data. 

The trends that we see in the XLFs of IC 348 directly re- 
flect the deficiency of X-ray luminosity in very low-mass disk- 
bearing stars responsible for the drop of the L^/ Lboi ratio to- 
wards fainter, less massive stars. We have compared this result 
to ONC samples constructed by combining the recent YSO clas- 
sification with COUP-derived parameters, and find a decline of 
the Lx/Lboi level towards low-luminosity for the diskless stars 
in the ONC. In disk-bearing stars of the ONC the scatter is large, 
possibly due to biases related to the difficulty of correcting the 
observed luminosities for high circumstellar extinction that is 
ubiquitous at the young age of Orion. 

Some scenarios for a decrease of Lx/Lboi level towards 
lower masses are: (i) transition from a solar-like to a convec- 
tive dynamo, combined with a lower efficiency of the latter one 
resulting in decreased X-ray production; (ii) smaller co-rotation 
radius for lower-m ass stars and the en suing centrifugal disrup- 
tion of the corona dJardine et al Mm) : (iii) a mass dependent 
fraction of stars with X-ray emission suppressed by accretion 
dFlaccomio et alJ2003l) . Concerning (ii), we note that indeed the 
IC 348 stars in our lowest mass bin have all short rotation pe- 
riods (< 5 d) while the period distribu tion is much broader fo r 
higher mass stars, an issue discussed bv lCieza & Balibe^ d2006l) . 
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Alexander et al. (in prep.) present a detailed study of the rota- 
tion/activity connection in IC 348 and find no significant differ- 
ences between stars in the saturated and in the super- saturated 
regime. Moreover, rotational effects should affect both disk- 
bearing and diskless stars in the same manner. Scenario (iii) is 
the only one that distinguishes between YSOs of given mass but 
in different evolutionary phases and would explain the different 
XLFs of very low-mass Class II and Class III objects. The search 
for a direct connection between mass accretion rate and X-ray 
luminosity could shed light on this question but, at present, it is 
hiimpered because of the small number of IC 348 members with 
accretion measures. In the available sample of 12 objects we find 
no relation between Mace and L^, except for the trend that can 
be ascribed to the well-known dependence of both parameters 
on mass. 

To summarize, we have examined one of the deepest avail- 
able X-ray data sets on a young stellar cluster. Thanks to its rich 
population of very low-mass cluster members IC 348 provides 
an excellent opportunity to study the X-ray activity level at the 
boundary and beyond the substellar limit. We have established 
for the first time with statistical significance the non-constant 
behavior of L^/L\-,oi for the very low-mass pre-main sequence. 
Although the sensitivity of the X-ray data is comparable to that 
of other recent studies of star forming regions such as XEST and 
COUP, the well-characterized cluster census way down into the 
BD regime results in a low X-ray detection rate for very low- 
mass objects and the low-luminosity end of the X-ray data of 
IC 348 continues to be dominated by insufficiently constrained 
upper limits. A much longer exposure would substantially in- 
crease the X-ray census at the low-mass end of the stellar se- 
quence in IC 348. Additional studies of accretion rates in IC 348 
would help to understand the interplay between accretion and 
X-ray emission. 
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Fig. 3 Merged Chandra/ AClS-l image created from the four observations listed in Table|2] IC 348 members in different YSO phases 
are color-coded as in Fig.[T] Cluster members without known YSO class are shown as black circles; many of these are secondary 
components in binaries. 
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Fig. 4 X-ray luminosity and median energy for the low-mass end of IC 348 versus spectral type. The color codes refer to different 
YSO classes (see Fig. [T]). The colored lines in the right panel represent the mean of the median energy for each YSO group 
iiTespective of spectral type and its standard deviation, i.e. including also higher-mass stars for which individual data points are 
outside the plotted range. 
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Fig. 5 Spitzer color-color diagrams for IC 348 members in the Chandra fields. X-ray detections are represented by filled circles 
and non-detections as x-points. The color codes for different SED YSO classes are the same as in the previous figures and, for 
convenience, defined agai n in the top left of eac h panel. The boxes in the left figure show the areas defining Class II according to 
iMegeath et alj ( |2004|) and lHartmann et al.1 (l2005h : dashed and sohd lines, respectively. 
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Fig. 7 X-ray versus bolometric luminosity for IC 348 (filled, thick plotting symbols) compared to the ONC (open, thin plotting 
symbols): (left) - Class n, (right) - Class III. Stars with indications for a disk warp provoking temporarily high extinction (AA Tau- 
Uke stars) are highhghted with large circles. 
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Fig. 9 Cumulative distributions of X-ray luminosity for different subgroups of YSO in three mass bins (from top to bottom). The left 
panels represent three YSO classes based on IR excess (red dotted - Class II, blue dashed - Class 11/111, and green solid - Class III) 
and the right panels represent YSO classes based on Ha emission (green solid - wTTS and red dotted - cTTS). 



15 



B. Stelzer et al.: X-ray view of IC 348 in the light of an updated cluster census 



1 — 1 — 1 — 1 — 1 — I — 1 — 1 — 1 — 1 — I — 1 — 1 — 1 — 1 — I — 1 — 1 — 1 — 1 — I — 1 — 1 — 1 — 1 — I j — 1 — 1 — 1 — 1 — I — 1 — 1 — 1 — 1 — I — 1 — 1 — 1 — 1 — I — 1 — 1 — 1 — 1 — I — 1 — 1 — 1 — 1 — I- 




log (L^/Lboi) log (L,/Lb„i) 




-5.0 -4.5 -4.0 -3.5 -3.0 -2.5 -2.0 -5.0 -4.5 -4.0 -3.5 -3.0 -2.5 -2.0 

log (L^/Lb J log (L,/Lb„i) 




Fig. 10 Same as Fig.|9]but for the fractional X-ray luminosity log (Lx/-^boi)- 
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Table 7. Primary Chandra Catalog: Basic Source Properties 





Source 




Position 










Extraction 










Characteristics 






Seq. No. 


CXOU J 


a (J2000.0) 


5 (J2000.0) 


Error 


e 


Ct t 


Ct t 


Bt 




PSF Frac. 


Signif. 




Anom. 


Var. 


Eff. Exp. 


-'-'median 






(°) 


(°) 


(") 


(') 


(counts) 


(counts) 


(counts) 


(cotints) 












fks) 


(keV) 




(2) 


(3) 


(4) 


(5) 


(6) 


(1) 






(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


fl6) 


(17) 


1 


034313.75+320045.2 


55.807303 


32.012563 


0.7 


10.4 


50.0 


9.5 


27.0 


27.5 


0.91 


5.0 


<-5 




b 


100.4 


2.1 


2 


034316.32+320305.4 


55.818019 


32.051512 


0.7 


10.8 


65.4 


10.4 


27.6 


54.4 


0.90 


6.0 


<-5 




a 


100.1 


3.5 


3 


0343iy.sy+320241.4 


55.832910 


32.044855 


0.8 


10.0 


27.9 


7.4 


17.1 


4.4 


0.90 


3.5 


<-5 






93.5 


1.3 


4 


034326.96+320042.1 


55.862341 


32.011705 


0.5 


7.8 


34.3 


7.1 


8.7 


3.5 


0.90 


4.5 


<-5 




a 


111.4 


1.4 


5 


034328.76+315458.9 


55.869862 


31.916371 


0.1 


7.5 


405.0 


21.0 


12.0 


241.3 


0.89 


18.9 


<-5 




c 


105.7 


2.5 


6 


034330.52+315649.6 


55.877202 


31.947123 


0.4 


6.6 


32.1 


6.6 


4.9 


25.6 


0.91 


4.5 


<-5 






94.3 


3.0 


7 


034332.08+320617.2 


55.883700 


32.104795 


0.4 


9.0 


111.7 


11.7 


11.3 


8.0 


0.89 


9.2 


<-5 




a 


67.3 


1.1 


8 


034336.33+320402.0 


55.901402 


32.067234 


0.6 


8.1 


27.6 


7.4 


17.4 


21.4 


0.91 


3.5 


<-5 




b 


172.7 


3.5 


9 


034336.47+320343.6 


55.901976 


32.062120 


0.5 


8.0 


47.2 


8.7 


16.8 


25.8 


0.89 


5.1 


<-5 




a 


173.4 


2.1 


10 


034342.05+315224.9 


55.925209 


31.873606 


0.5 


6.9 


19.5 


5.7 


6.5 


11.8 


0.90 


3.1 


<-5 




a 


111.9 


2.7 



in 



Note. — Table|2]is available in its entirety in the electronic edition of A&A. Interesting sources mentioned in the text are shown here for guidance regarding its form and content. ^. 

Col. (1): X-ray catalog sequence number, sorted by RA. ^ 

Col. (2): lAU designation. o 

Cols. (3) and (4): Right ascension and declination (in decimal degrees) for epoch J2000.0. hh 

Col. (5): Estimated standard deviation of the random component of the position error, y'trl+o^. The single-axis position errors, Ux and cry, are estimated from the single-axis standard deviations of the PSF inside thei>j 
extraction region and the number of counts extracted. 

Col. (6): Off-axis angle. 5' 

Cols. (7) and (8): Net counts extracted in the total energy band (0.5-8 keV); average of the upper and lower la en'ors on col. (7). & 
Col. (9): Background counts expected in the source extraction region (total band). 
Col. (10): Net counts extracted in the hard energy band (2-8 keV). 

Col. (1 1): Fraction of the PSF (at 1.497 keV) enclosed within the extraction region. A reduced PSF fraction (significantly below 90%) may indicate that the source is in a crowded region. q 
Col. (12): Photometric significance computed as net counts divided by the upper error on net counts. 

Col. (13): Logarithmic probability that extracted counts (total band) are solely from background. Some sources have Pg values above the 1% threshold that defines the catalog because local background estimates can rise^ 
during the final extraction iteration after sources are removed from the catalog. 

Col. (14): Source anomalies: (g) fractional time that source was on a detector (FRACEXPO from inkarj) is <0.9. 
Col. (15): Variability characterization based on K-S statistic (total band): (a) no evidence for variability (0.05 < Pks)'^ (b) possibly variable (0.005 < Pks < 0.05); (c) definitely variable (Pks < 0.005). No value isp, 

reported for sources with fewer than four counts or for sources in chip gaps or on field edges. 2- 

Col. (16): Effective exposure time: approximate time the source would have to be observed at the aimpoint of the ACIS-1 detector in Cycle ??? to obtain the reported number of net counts (see ... S 

Col. (17): Background-corrected median photon energy (total band). ^ 



